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Abstract 

A class of second order difference (discrete) operators witli a par- 
tial algebraization of the spectrum is introduced. The eigenfuncions 
of the algebraized part of the spectrum are polinomials (discrete poli- 
nomials). Such difference operators can be constructed by means of 
Uq{sl2), the quantum deformation of the SI2 algebra. The roots of 
polinomials determine the spectrum and obey the Bethe Ansatz equa- 
tions. A particular case of difference equations for g-hyper geometric 
and Askey- Wilson polinomials is discussed. Applications to the prob- 
lem of Bloch electrons in magnetic field are outlined. 
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1 Introduction 



In this paper we attempt to describe difference equations of second order in 
one variable, 

a{z)^{qh) + d{z)'^{q~h) + v{z)^{z) = E^{z), (1) 
having polynomial eigenf unctions of the form 

N 

<iliz)=l[iz-zj. (2) 

m=l 

Here g is a parameter and the polinomials are parametrized by their roots 
Zm- We call it algebraization of the spectrum. 

Difference equations are closely related to the discrete equations. The 
discrete equation 

ttni^n+l + dntpn-l + Vnlpn = Elpn ■ (3) 

may be obtained from the difference equation (|1|) by setting z = g^": 

vl>(g2«) = a(g2«) = a„, rf(g2«) = 4, ^;(g2n) = . (4) 

Therefore, if the difference equation has a polynomial solution, then one can 
find a solution of the corresponding discrete equation as a polynomial on a 
discrete support: 

N 
m=l 

Generally only a part of the spectrum of difference and discrete operators 
can be algebraized. However, in some cases the entire spectrum is algebraic. 
The important examples of this type are periodic discrete equations. They 
appear if g is a root of unity, i.e., g^*^ = 1, where Q is an integer. Then we 
obtain equations with periodic coefficients a„ = an+q, dn = dn+q, Vn = Vn+q 
and the periodicity condition 

t/jn = i)n+q- (6) 

Some of the equations (|l],^ have important physical applications. In par- 
ticular, the method of algebraization has been applied [jl], to the problem of 
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Bloch electrons in magnetic field on a lattice (the Azbel-Hofstadter problem) 

Discrete equations with periodic coefficients are reacher than the differ- 
ence equation with |g| = 1. One can impose the quasiperiodic boundary 
condition 

= e^'>n+Q, (7) 

where k is the Bloch momentum. In this case the spectrum generally has 
Q bands. The periodic boundary condition describes the crossection of 
the band spectrum at = 0. The quasiperiodic equations {k ^ 0) have no 
difference analog. We do not consider them in this paper. 

The class of algebraized operators may be extended by the " gauge" trans- 
formation '^{z) U{z)^>{z), where U{z) may not be a polynomial. Then 
a{z) U{q^z)-^a{z)U{z), d{z) U {q-^ z)-^ d{z)U (z) , v{z) v{z). In par- 
ticular, one can always choose the Jacobi operator to be symmetric: 
a{z) = d{q^z). 

In Appendix A we show that the difference equation (|I]) has polynomial 
solutions if a{z)^ d{z), v{z) are certain Laurent polinomials of order 2. They 
are determined by 7 parameters (not including q) and may be computed di- 
rectly. However, the algebraization of difference equations has an intimate 
relation with the representation theory. The relation between representations 
of the SI2 Lie algebra and second order differential equations having polyno- 
mial solutions is known ( see e.g.[^ ^, ^, for early works, ^, |TD[ for recent 



systematic treatment and ||Tl| for a review). In this paper we show that dif- 
ference equations having polynomial solutions can be classified according to 
representations of the quantum algebra Ug{sl2), the g-deformation of 5/2- 

To construct linear operators with partially algebraized spectrum we em- 
ploy the following strategy. Consider the algebra V of difference operators 
acting in the space of complex functions. Let A C V he its subalgebra such 
that A has a finite-dimensional irreducible representation. This means that 
elements of A leave invariant a finite-dimensional functional subspace. Some 
of the eigenfunctions of operators from A must then belong to this subspace. 
Let us choose the invariant functional subspace to be the linear space Poln 
of polinomials of degree at most n. Our task, then, will be to realize a finite 
dimensional representation of the algebra A in this space. 

The basis of the space Poln may be choosen as monomials z^ , k = 
0, 1, 2, . . . , n. In this basis, the algebra A is naturally decomposed into the 
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raising (lowering) parts {A-) and diagonal operators (Cartan subal- 
gebra) with the property A-^-Aq C A+, A-Aq C A-. The space Poln is in- 
variant if and only if A+z^' = A-Z^ = 0. This means that the representation 
of A must have highest and lowest weights. Then diagonalizable operators 
from A must have n + 1 different polynomial eigenfunctions. When the sub- 
algebra A and its representation are given we want to select the difference 
operators of second order. In this case roots of the polynomial eigenfunctions 
satisfy Bethe equations. 

This strategy has been used for construction and classification of differ- 
ential operators having polynomial eigenfunctions 0, ^ |, |T0|, 0, |12 



different approach was suggested in 0). Let us use here the same notation 
V for the algebra of differential operators. Then the subalgebra A G V is a 
factor of U{sl2) (the universal enveloping of s/2) over its center. The algebra 
SI2 can be realized by first order differential operators, 

In the invariant subspace Pol2j the representation (H) has highest and lowest 
weights ±j (integer or half-integer spin of the representation). This realiza- 
tion provides the embedding of A into the algebra of differential operators 
T>. More precisely,^ in this approach is identified with the factoralgebra 
U{sl2)/{S'^ — j{j + 1)) over the ideal generated by S"^ — j{j + 1), where S"^ is 
the Casimir (central) element. 

Then the Hamiltonian of the Euler top, 

H = J2a,,S,S,+J2(3,S,, (9) 

i.e., a general bilinear form in the s/2 generators gives a family of second 
order differential operators having in general 2j + 1 independent polynomial 
eigenfunctions. 

Similar arguments may be applied to the quantum algebra Uq{sl2)- In 
this case one may identify A with a certain factor of Uq{sl2) over its center. 
Its highest and lowest weight representations ( g-analogs of (§)) are known 
and presented in Sect. 2. In this paper we show that homogeneous bilinear 
and linear forms in generators of Uq{sl2) (g-deformations of the Euler top 
(ID) provide the difference equations of interest (Sect.3). 
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In addition to Uq{sl2) some other deformations of the s/2 algebra are 
known (Appendix B). They also can be used to construct difference equa- 
tions with polynomial solutions (we note that a deformation not equivalent 
to Uq{sl2), which generates a subalgebra of Ug{sl2), has been used in the 
paper in attempt to construct algebraized difference operators (see Ap- 
pendix B)). Among different deformations the quantum algebra Uq{sl2) plays 
a special role. It indicates a relation of difference algebraized operators and 
classical discrete integrable systems (in fact the former are Lax operators 
for a class of nonlinear integrable equations). In this paper we do not use 
the coproduct structure of the quantum algebra. However, we anticipate an 
upgrading of the present approach for quantum integrable systems. 

Recently we have shown that the spectrum of the Harper equation (also 
known as the discrete Mathieu or "almost Mathieu" equation) 

iJn+i + + 2 cos(A; + n^)ip^ = Eip^ (10) 

at cos{kQ) = — 1, where Q is the denominator of the parameter $ = 27rP/Q, 
is algebraized 0: 

m=l 

The roots obey the Bethe Ansatz like equations 
and determine the spectrum []: 

m=l 

The Harper equation describes the Bloch particle on a square lattice in mag- 
netic field (Azbel-Hofstadter problem). In this paper we write the Bethe 
Ansatz equations for a general class of difference operators having polyno- 
mial eigenfunctions. In addition, we present a solution for the Bloch particle 
on the triangular lattice in Sect. 4. 

more general algebraization in terms of meromorphic functions on algebraic curves 
rather than polinomials has been introduced by Faddeev and Kashaev for quasiperiodic 
Harper's equation for an arbitrary k. 
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In Sect. 2 we give main formulas related to Uq{sl2) and its representations. 
In Sect. 3 we consider the spectral problems for difference operators appear- 
ing as linear and bilinear forms in generators of Uq{sl2)- They are difference 
(discrete) analogs of the algebraic forms of known differential equations hav- 
ing polynomial solutions ||, ^ |10|, |Tl|, [l^, ^ . Besides, we point out a class of 



difference spectral problems solvable in terms of g-deformed classical orthog- 
onal polinomials with discrete measure (big g-Jacobi polinomials). Possible 
applications to the problem of Bloch particle in magnetic field are outlined in 
Sect. 4. In Sect. 5 the continuum ("classical") limit is discussed. A related 
class of difference equations with solutions of the form of symmetric Laurent 
polinomials is treated in Sect. 6. The solutions of a particular case of these 
equations are Askey- Wilson polinomials. 



2 The quantum algebra Uq{sl2) and difference 
operators 

The algebra Ug{sl2) (a g-deformation of the universal enveloping of the s/2) 
is generated by the elements A, B,C, D, with the commutation relations [|16|, 
^, 

AB = qBA, BD = qDB, 
DC = qCD, CA = qAC, 

AD = 1, [B,C] = -. (12) 

q — q^^ 

The deformation parameter q may be considered as a formal variable. In the 
classical limit q —>■ 1, the quantum algebra turns into the universal enveloping 
of sh: {A - D)/{q - q-') -> S3, 5 ^ S+, C ^ 

The central element of this algebra is a g-analog of the Casimir operator 

n^t^^^^BC. (13) 

(g-g 

As g 1, n - 2(g - q-^y^ tends to S"^ + 1/4. 
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The commutation relations (|T2|) are simply another way to write the 
intertwining relation for the L-operator: 

R{u/v){L{u) ® 1)(1 ® L{v)) = (1 ® L{v)){L{u) ® l)R{u/v) (14) 

with the trigonometric i?— matrix 

R{u) 



1 



1 



\ ^) + -{q - l){u + q \ ^)a3 
+ {q - q'^){a+ (g) a-_ + (T_ (g) cr+) 

satisfying the Yang-Baxter relation {aj are Pauli matrices; a± = 
Generators A, B, C, D are matrix elements of the L-operator 

L{u) 



q-q- 

B 



C 

ukD-u-^k-^A 



q-q- 



5 + 

(15) 

:(Ti±za2)/2). 
(16) 



Here u is the spectral parameter and k is an additional parameter ("rapidity" 
at the site). Note that the i?-matrix is the L-operator in the spin 1/2 rep- 
resentation. It is given by the same matrix (|I^) for k = q^l"^ with elements: 
A = q^''^, D = q-^"'-', B = a+, C = a^. 

Irreducible finite dimensional representations of dimension 2j + 1 can 
be expressed in the weight basis, where A and D are diagonal matrices: 
A = diag {q\ q~^). An integer or halfinteger j is spin of the representation. 
The value of the Casimir operator ([1^) in this representation is 



(q-q 1)2 



(17) 



These representations can be realized by difference operators acting 
in the linear space of polynomials F{z) of degree 2j. Let us introduce "shift" 
operators and T_: T±F{z) = F[q^^z), T^T_ = 1, 



q^h-'T. 



(18) 



Then 



A = q~^T+, D = q^T^, 

B = z{q-q-^)-^ (q^^T_ - q-^^T^ 



C 



(T_ - 



(19) 
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Then the lowest weight vector is Fq{z) = 1 whereas F2j{z) = z^^ is the 
highest weight vector: CFq = , BF2j = 0. The reahzation ([T9| ) is a smooth 
g-deformation of the representation of the s/2 algebra by first order differential 
operators (||). The Casimir operator (|T3|) in this realization is equal to the 
c- number Qj ([T7|) . This means that ([l9| ) actually gives a representation of the 
factoralgebra Uq{sl2)/{fi — fij) over the two-sided ideal generated hj fl — flj. 

If g is a root of unity there is, in addition, three parametric family of 
finite dimensional representations having, in general, no lowest and no highest 
weight ||2^ . Sometimes they are called cyclic (or unrestricted) representations 
np. They correspond to discrete quasiperiodic equations (0) which has no 



direct difference analog 



3 Second order difference operators related 

to Uq^sk) 

We call a difference operator in a variable z any linear combination of the form 
J2i fi{,z)Tl^, where fi{z) are rational functions of z and the sum is finite. Such 
operators form the algebra V mentioned in the Introduction. The general 
form of a difference operator of second order then is (cf. (|I])) 

A = h{z)Tl + Uz) + h{z)T^^ , 

where s is an integer. 

Second order difference operators appear as linear (s = 1) or bilinear(s = 
2) forms in Uq{sl2) generators. In fact operators coming from linear forms are 
equivalent to some bilinear forms due to isomorphism of the algebras gener- 
ated by A, q^/'^BA, q^^'^CD, D and A, B, C, D. Nevertheless, it is convenient 
to consider them separately. We begin with linear forms. 

3.1 Linear forms 

Consider a linear form in the quantum algebra generators 

L = aA + dD + {q- q-^) {bB + cC) , (20) 
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where a, b, c, d are parameters. The diagonahzation of this operator in 
2j + 1-dimensional representations of Uqi^sh) leads to the difference equation 

{cz-^ + aq-^ - hq-^^z)m{qz) + {-cz'^ + dq^ + hq^h)m{q-^z) = E^{z). (21) 

This equation has 2j + 1 polynomial solutions. Let us parametrize a poly- 
nomial by its roots (0). Pluging (||) in (^) and dividing both sides by "^{z) 
we get 

N 2 N -2 

n '-^^ + ^(-) n = E (22) 

where we have denoted 

a{z) = cz-^ + aq-^ - bq-^h , (23) 

d{z) = -cz-^ + dq^ + bq^^z. (24) 

Assume that c 7^ 0. Then two different cases are possible: (1) 6 7^ 0, c 7^ 
0, (2) 6 = 0, c 7^ (the case 6 = c = is trivial). 

Let us first consider the general case (1). Assume that '^{z) is nondegen- 
erate, i.e., all different. The l.h.s. of (0) is a meromorphic function, 

whereas the r.h.s. is a constant. To make them equal we must cancel all the 
singularities of the l.h.s. They appear at singular points of a{z), and d{z) 
(simple poles at z = and z = 00) and ai z = z^- The singular part at 2; = 
vanishes automatically (note that no one of be equal to 0). The 

residue at infinity is equal to bq'^^~^{l — g^^"^-'). Its vanishing determines 
the degree of the polynomial: N = 2j. Comparing the constant terms in the 
both sides of ( p2|) one finds the energy spectrum: 

2i 

E = aq^ + dq-^ - b{q - q-^) J2 Zm- (25) 

m=l 

Annihilation of poles aX z = z^ gives the following Bethe-ansatz equations: 

bzf + dq-^zi - cg~^^ ^ _ ^ qzi - z^ , . 

bq-'^hf - aq-^zi - c ^ Zi - qz^ ' 

In the case (2) the operator L is triangular in the basis of monomials. 
It includes only A,D and C (the generators of the Borel subalgebra). As a 
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result the coefficients of "^{z) can be recursively determined. Therefore there 
is only one polynomial solution for each degree N = l,...,2j. One can see 
it from (P^). Its l.h.s. is regular at z = oo hence there is no restriction on 
degree of the polynomials: 

dq^i - c ^ SL qzi- Zm 

^ 11 Z —■ (27) 



aq ^Zl + C m=l,m^l - l^r. 

The energy spectrum is given by the simple formula 

= aq^-^ + dq^-^. (28) 

The Bethe equations (|27|) must have exactly one solution (modulo permuta- 
tion of the roots). In section 3.2 we discuss a relation of the solutions with 
g-orthogonal polynomials. 

3.2 Bilinear forms 

Consider a general bilinear form in Uq{sl2) generators: 

G = aA^ + dD'^ + (g - q-^){c2CA + b2BD + b^BA + c^CD) + 

Hq-q'')\hB^ + ciC^). (29) 

where a, d, Ci, hi {i = 1,2,3) are arbitrary parameters. 

In the representation of spin j it is a difference operator: 

G^{z) = a{z)^{q'^z) + d{z)'^{q-'^z) - v{z)^{z), (30) 

where 

a{z) = hq~^^+^z^ - hq-^h + aq-^^ + C2q'^ z'^ + c^q'^ z-^ , (31) 

d{z) = hiq^^-^z^ + h2q^^z + dq^^ - c^q^ z-^ + Ciqz-^ , (32) 
v{z) = {q + q-'){b,z^ + ciz-^) + (csg'^' - C3g^>-^ + (fesg"^' - hq')z . (33) 



The difference operator given by the linear form ( pi|) is a particular case of 
the bilinear form 

aA^ + dD"^ + {q- q-^){q^cCA + q^bBD + q-^bBA + q-^cCD) , 
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where q is to be changed to q^. 

Let us plug d^) in (|30D and divide both sides by '^{z). We get 

N 2 N —2 

a(z) n ^ '~ + d{z) n ^-4^ - viz) = E . (34) 

Z Zm Z Zfn 

Suppose that at least one of the coefficients ci, C2, C3 is nonzero, then 
there are three main different cases: 

(i) at least one of bi, 62, &3 is nonzero and both a{z) and d{z) (see (|3T|) , 
(P^) are nonzero (this is the case of generic position; the other two can be 
considered as exceptional cases); 



(ii) All 6's are zero: the quadratic form ( pQ] ) includes only A, D and C 
(generators of the Borel subalgebra of Uq{sl2)). In this case the eigenfunc- 
tions are big g-Jacobi polynomials [^. They include all g-deformed classical 



orthogonal polynomials with discrete measure |2l 

(iii) One of the functions a{z) or d{z) is identically zero. In this case the 
second-order difference equation ([30| ) reduces to a ffist-order equation. 

Let us consider the case (i) ffist. The l.h.s. of (0) is a meromorphic 
function, whereas the r.h.s. is a constant. To make them equal we must 
cancel all the singularities of the l.h.s. They appear at singular points of 
a(z), d{z) and v{z) (double and simple poles at z = and z = 00) and at 
z = Zm (again, we consider non-degenerate case when all of them are simple 
poles). Note that the case when at least one of Zm^s is zero (i.e., coincides 
with the pole of a{z) and d{z)) needs special consideration. At the moment 
we assume that Zm ^ 0- The singular part a.t z = vanishes automatically. 
Vanishing of the singular part at 2; = 00 , 

N 

bsiq' - g'^-'^> + z{q - q-'Mq'""-'^ - q-'""^'^) z^ , (35) 

m=l 

determines degree of the polynomial: N = 2j. If 61 = and 62/^3 = —q^^ 
for an integer M > —j , N = j + M is also possible. Below we consider only 
the generic case N = 2j. 

Comparing the constant terms in the both sides of ( ^4]) we find the energy 
spectrum: 

E = bi{q - q~^){q'^ - q'"^) Y ZnZm - 
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2i 

~{q - q-'){b2q-'^' + hq'-') £ + + dq'"' • (36) 

m=l 

Finally, annihilation of poles ai z = Zm gives the following Bethe-Ansatz 
equations 

^-q'^ h /-I,. ..,2,. (37) 

The Bethe equations is a system of 2j algebraic equations. It must have 
exactly 2j + 1 solutions corresponding to different eigenfunctions. In the 
case (i) all of them are polynomials of one and the same degree 2j. 

Note that d{z)/a{z) is a rational function having at most 4 zeros and 4 
poles. Therefore ( P7| ) looks like a system of Bethe equations for a XXZ-spin 
chain on at most 4 sites with different spins at the sites. The reader familiar 
with the algebraic Bethe ansatz should notice that eq. is the 



Baxter identity for eigenvalues of the transfer matrix t{z) for this system 

^{z)t{z) = a{z)'^{q'^z) + d^z)"^ (q'^ z) , (38) 

In this context the eq. ([38|) determines t{z) and "^{z) provided that a{z) and 
d{z) are known functions. 

In the special case q'^^'^^ = 1 the difference operator has finite dimension 
2j ' + 1. Therefore polynomial egenfunctions and the Bethe Ansatz cover all 
the spectrum. There are no more solutions other than polynomials. 



3.3 Triangular and first order operators 
3.3.1 g-Hypergeometric equations 

The difference operators that include only elements A, D and C (case (ii)) 
lead to g-hypergeometric equations. These operators preserve not only the 
space Pol2j for some j, but all the spaces Po/„ for any n > 0. Indeed, they 
are lower triangular in the basis of monomials z'^. The l.h.s. of ( p4|) is now 
regular at 2; = oo , so there is no restriction on degree of the polynomials. 
The Bethe equations are valid for any < 2j ' + 1: 

^ = q^^ n / = !,..., AT. (39) 
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Triangularity leads to a very simple structure of the spectrum of G\ 



E 



N = aq 



2N-2j 



+ dq'^-'\ iV = 0,...,2j. 



(40) 



A general operator of this kind has the same spectrum as the diagonal oper- 
ator aA^ + dD^\ it does not depend on the coefficients ci, C2, C3. 

For the sake of completeness let us identify the solutions of (|30|) in the 
case (ii) with known g-orthogonal polynomials. Consider the big g-Jacobi 
polynomials P^'^\z]'^,5]q^) (where we have used the standard notation 
from the textbook ^3[; to shorter the notation we denote them as P„(z)). 
They can be expressed in terms g-hypergeometric series: 



P^'3)(z;7,5;g') =3 02 



^-2n „2Q+2/3+2n+2 „2a+2^,-l^ 



(41) 



and obey the difference equation 

+ (1 + {5- 7)^-^ - ^5z-^)Pn{q-^z)- 
-[((1 + g2")5 - (1 + _ ^5(1 + q~^)z-^]P^{z) 

= (g-2" + g2n+2a+2/3+2)p^(^) _ 



(42) 



This equation is equivalent to (|3 

6 = C2,d~^q~^ , —75 = cid 



7 - 

ad~ 



provided d ^ 0, bi 

2a^ 



and 



• 1 



g2/3^ 



1^-2- 



4i 



Note that (^) is formally symmetric with respect to the change n 
—a — (3 — n — 1. This symmetry becomes important ifm = a + /5 + n + lisa 
negative integer. Then the series (^) truncates at ra-th term if n < \m\, but 
if |m| < it truncates at |m|-th term. One has a + jS = \ogq2{a/d) — 2j — 1. 
In particular, if a = ci the series (^Tj) is symmetric with respect ton ^ 2j — n. 
In other words, for integer j it gives only j + 1 different eigenf unctions (they 
are polynomials of degrees 0, 1, j). It can be easily seen that in this case 
the original operator contains Jordan cells and, therefore, is not completely 
diagonalizable (i.e., it has less than 2j + 1 different eigenfunctions) . In this 
case the formula (^Tj) still gives all the eigenfunctions. 

At (5 = 0, the big g-Jacobi polynomials Plf'^\z;'y,0;q) reduce (up to 
an ra-independent factor) to the little g-Jacobi polynomials ^^^'"^(^7"^, g). 
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They provide eigenf unctions of the hnear form 



at 6 = 0. In this case 



qh 



-ad' 



-iq-i-2j 



7 



cd K It is known that the little q 



Jacobi polynomials can be expressed through the g-hypergeometric function 



2V1- 



pi^^"\zr\q) =2 01 



H 1 H 



q 



/3+1 



(43) 



Again, if the operator is not completely diagonalizable, ( ^5] ) gives less than 
2j + 1 different eigenf unctions. 

The following comment is in order. In contrast to the linear forms (^) the 
system of algebraic equations (|39[) generally has more than one solution for 
a given N . However, only one of them corresponds to an eigenfunction. The 
point is that the Bethe equations have been derived under assumption that 
all roots are nonzero. In fact one or more zero roots are not forbidden (for 
instance, g-Legendre polynomials of odd degree are odd functions). However, 
they change the residues at 2r = of the first two terms in (@), so this case 
needs a special consideration. To illustrate this let us give an example. 

Example. Consider the spectral problem for the following quadratic forms: 



Gi = A' - + {q- q-'){q^-'CA + BD-BA- q'-^CD), 

G2 = q-'A^ + qW^ + (g - q-^){q-'CA + q^CD) + (g - q-^fC\ 

Though Gi corresponds to the case (i) and G2 to the case (ii) they lead to 
identical Bethe equations. For instance, at j = 1/2 the equation is 

qzl — g^/^Zi + g^/^ 
q-^zl + q-y^zi + g-V2 = 

It has two roots: = and , z'l = 2(g^/^ — g^^/^)^^. Both of them give an 
eigenfunction of Gi, while for G2 z[ is an artifact.) 



3.3.2 First order operators 

At last, we consider the case (iii).| When, say, d{z) = the operator is 

Go = aA^ + {q-q-'){c2C + hB)A. (44) 
^ These operators appear in harmonic analysis on the quantum group |26[, p7|. 
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The Bethe equations (^) are simplified as 



2j 

«(^«) n (^^^' - =0, 

k=l,k^l 

SO eigenf unctions and the spectrum can be easily found. Let z± be the 
two roots of the quadratic equation a{z) = 0, then the m-th eigenfunction 
(m = 0, 1, 2j) and the spectrum are 

2j — 1— m m—l 

^m{^)= n - n - , 

k=0 1=0 

Em = hq-^iz^q^^-^"' + ^+g2--2i) _ 

3.3.3 Algebraic structure unifying triangular and first order op- 
erators 

The similarity of the spectra of operators of types (ii) and (iii) is not an 
accident. It turns out that these operators form a simple quadratic algebra. 
Consider two operators of the type (iii) 

Hi = fxiBA + fi2CA + /igA^ , (45) 

H2 = viDB + V2DC + z/gD^ (46) 
with arbitrary coefficients and z/j and a general operator of the type (ii) 

Hz = [q -q )l _ g/i2Z^2C - /i2Z^3-DC - /i3Z/2C A I . 

^ ^ (47) 
It is straightforward to check that these operators form the closed algebra 

q-^H,Hj - qHjHi = g^Hk + hu (48) 

where {ijk} stands for any cyclic permutation of {123}. The structure con- 
stants Qi and hi may be expressed in terms of coefficients /ij and Vi and the 
value of the Casimir operator Vt (|13|). 

A similar algebra has been considered by Granovskii and Zhedanov in [ p8[ , a particular 
case hi = h2 = = has been discussed in pOt], p2[. 
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The algebra ( ^S]) allows one to find the spectrum of i^i, H2 or in an 
algebraic way Suppose one knows an eigenvector of Hi. Then, some 

linear combination of Hi is a creation operator: it creates a new eigenvector 
by acting to the known one. As a result the spectrum of Hi is found to be a 
trigonometric function of the number of level. 



4 Difference periodic equations and the group 
of magnetic translations 

Difference periodic equations and discrete equations with periodic coefficients 
appear when g is a root of unity: 

q = exp(mP/g), (49) 

where P and Q are coprime integers. Some of them have important applica- 
tions in physics. 

Let us briefly recall the problem of Bloch particles in magnetic field 
(Azbel-Hofstadter problem). Consider a particle on a two-dimensional square 
lattice. The Schrodinger equation has the form 

J2 t;^e*^'^("V(^ + /^) = E^Pin) , (50) 

where n is a lattice site, /i is a lattice vector and tp = is a hopping am- 
plitude. The gauge potential A/j(n) = —A_pi{n — p.) describes a homogenious 
magnetic field with flux $ = l-aPlQ per plaquette of the square lattice: 

plaquette 

In this context the parameter q has a clear interpretation: it is a flux per 
half of the plaquette. The problem may be expressed in terms of magnetic 
translations ||^: 
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Let X and y be unit vectors along x and y directions on the square lattice 
and n = {rix^ny). They generate the algebra of magnetic translations 

'~r^'~r^ ^nym^—nxniyrp^ ^ rp—l ^ 

-'-n-'-m y -'n+m; ^ J- —n i 

TyT^ = q^T^Ty, TyT_^ = q~^T_^Ty. (51) 

The operators T? are central elements. As an algebra over its center the 
algebra of magnetic translations is the finite Heisenberg-Weyl algebra. 
In these terms the Hamiltonian of the problem is 

H = T.tf^T^- (52) 

The quantum algebra Uq{sl2) (|12|) may be expressed through the Heisenberg- 
Weyl generators, i.e., in terms of magnetic translations There are many 
different ways to do it (see e.g. [BH] ). One of them (for odd P) is 

BD = iq^'^q - q-^)-\T^ + T_y), CD = tq-^l\q - q-^)-\T.^ + T_2x-y), 

BA = iq-^l\q - q-')-\T^ + T^^+y), 

B' = -iq- q-')-\T^-y + Tg.+y + (y + y-^)r2x), 

C^ = -{q- q-')-HT-.+y + T_3x-y + (y + y-^)T_2x). 

(53) 

The central element (|T3p in this representation is 

n = -2(y - y-i)-2. (54) 



Comparing with ( |T7| ) we see that this value corresponds to the representation 
(pJl) with j = (Q — l)/2 and dimension Q. The same value (0) corresponds 
to a 2-parametric family of the Q -dimensional cyclic representations . 
Pluging (pi) into (pOT) we obtain 



G — CtTx+y + (iT_x-y — ^l^x-y — ClT_x+y — &1^3x+y — ClT_3x-y + 

+2yi/2(c2Ty + 62r_y) + «yV2(62rx + C2r_x) + tq-'^HhT^ + C3r_x) + 

+Zy-l/2(63T2x+y + C3T_2x-y) " (? + g"')(ClT_2x + &lT2x). 

(55) 
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The operator G is hermitian ii d = a, hi = ci, —qh2 = C2, &3 = —qc^- It 
describes a particle in a magnetic field on the square lattice with a hopping 
along X — y and x + y diagonals, along x, y and 2x bonds and along 2x + y 
and 3x + y diagonals. 

Although the Hamiltonians (|55D contains many parameters, the number 
of interesting problems is limited since the hopping amplitudes must not 
depend on the magnetic flux (i.e., on q). Then only few cases remain. Among 
them, there is a triangular lattice: 61 = Ci = 63 = C3 = 0, C2 = &2 = i^g"^^^, 
a = d = tj^O (the sign depends on parity of (P-l)/2: it is "+" for (P-l)/2 
even and "— " otherwise), 

H = t{T^ + y + T_,_y) + + + Ty+ T„y . (56) 

Here the triangular lattice is treated as a square lattice with x + y diagonals. 
The square lattice may be obtained at t = 0. This case has been considered 
ini. 

In terms of generators of Uq{sl2), the Hamiltonian ( p6D is 

H = t{A^ + D^) ± tq-^/\q - q-'){CA + BD) (57) 
and under the representation (|T^) we obtain the spectral problem 

{z-^ - tq)'^{q'^z) + (zq-^ - tq-^)"^ (q-"^ z) - {z + z'^)'^{z) = E^(^) (58) 

How does this equations appears from the original problem? Let us choose 
a specific gauge to implement the flux $/2 per elementary triangle: 

A, = -$n^, Ay = $n„ A,+y = $/2. (59) 

Then, the Schrodinger equation aquires the form 

+e'*"-V^(n„ ny + l) + e-'^'^V^K, Uy - 1) (60) 
+t(e^*/2V;(n, + 1, + 1) + e'^'/^^in, - 1, - 1)) = E^{n,, riy) . 

Solutions are the Bloch waves: 
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where k = {k^, ky) is the Bloch momentum and ipn is Q-periodic discrete 
function. In these terms (^OD turns into the following version of the Harper 
equation: 

^^ikx-i^n _j_ ^^i{kx+ky)+i^/2^^^^^ _j_ ^^-ikx+i^n-i^ _j_ ^^-i(kx+ky)-i^/2'^^^ _^_|_ 

+2 cos(n<l> + ky)lpn = Eipn- 

(61) 

The spectrum of the problem has Q bands E = Ei{kx, ky), i = 1, Q. For 
kx = 0, ky = n ("midband points" of the spectrum) Eq.(pl|) becomes 

_ tq)i;^^, + (g2-2 _ tq'^)^^., - (g^" + g-2-)^„ = E^,, (62) 

with the periodic boundary condition @. This discrete equation is equivalent 
to the difference equation ([58|). Applying the results of Sect. 3 to (0)'® 
we obtain the Bethe equations 

Zi{zi-tq) _ 9^ q^zi - Zm , . 

+ m=l,jtl Zl-q Zm 

and the energies of the mid points of each band: 

Q-i 

E={q^-l)Y. Zm-t{q + q-'). (64) 

m=l 

Harper's equation (^) for an arbitrary Bloch momentum k may be ob- 



tained from { p7\j by using cyclic representations These representations 
generally has no highest or lowest wieght. The wave functions are not poly- 
nomials in this case. Nevertheless the Bethe Ansatz a kind of algebraization 
in terms of meromorphic functions on higher genus algebraic curves is still 
possible. The method which generalizes the algebraic Bethe Ansatz has been 



developed in Ref.[^ and applied to the chiral Potts model. Recently, Fad- 



deev and Kashaev |T^ used this method for the Harper equation on the 



square lattice to obtain the Bethe Ansatz equations for arbitrary k. 

Simplicity of the spectrum (pO]) of triangular operators considered in the 
Sect. 3 has a clear interpretation in terms of magnetic translations. In this 
case the Hamiltonian allows hopping only to the South- West, North-East and 
North- West directions on the lattice. Therefore, no closed loop trajectory is 
possible. As a result, the particle does not feel the magnetic field and the 
spectrum remains free. 
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5 Continuum limit. Differential equations 



Difference equations of Sect. 3 may liave different continuum limits. Let us 
consider the simplest one, when we assume that {<iJ{q^z)-^{q~^z))/{q-q-^) 
has a regular limit as g 1. Let us set q = and also assume that 
q^ = 1 + 0{h), so j is a fixed integer or half-integer number (let us note that 
this condition is not valid for the periodic difference equations considered in 
the Sect. 4, where q'^^'^^ = ±1). Then 

A = i + hS3 + h'^s^^ + oih^), D = i-hS3 + h^sl + o{h^), 
B = s+ + o{h^), c = S- + o{h^) 

as ^ — ^ 0. Here S±, are generators of U{sl2) realized as differential 
operators (|]). 

The coefficients of G (|29|) may also depend on Ti. Suppose they are regular 
at = 0. Up to the first order in % one has: 

a = a'^^^ + ha^^\ d = +hS^\ h, = hfHhhf\ Ci = cfHhcf\ z = l,2,3. 

In order to get a non-trivial limit one should put a^^^ = = —cf\ 

62"'* = —bf^- Then we obtain a general bilinear form in the s/2 generators 
(the quantum Euler top (P)). Components of the matrix aij and the vector /3i 

are 033 = a^^\ a++ = 2bf\ a = 2ci'\ a_^3 = 2c2 \ a;+_3 = — 262°\ «+- = 

a3,+ = as,- = 0, /33 = a^'^ - (3. = c^^ + c^^\ (3+ = h^^ + . Clearly, a 
global 5'L(2)-rotation leaves the spectrum invariant. Therefore, in the case of 
generic position, one may diagonalize the matrix aij and reduce the number 
of parameters to 4. Indeed, there are 3 momenta of inertia (eigenvalues of 
aij) and direction of the "magnetic field" (3 (parametrized e.g. by 2 Euler 
angles). However, the Casimir element S*^ contributes only to a c-number 
term, so we are left with 4 parameters. In terms of differential operators the 
adjoint action of the group amounts to linear fractional transformation. We 
do not know g-analog of this transformation in the discrete case. An explicit 
form of the differential equation which corresponds to the quantum Euler top 
in the representation (||) is 

Q4(^)^ + {Q2{Z) - (j - + (IjU - l)Ql{^) - JQ2{Z)) , (65) 
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where Qk{z) are general polinomials in z of degree at most k. These are the 
only differential operators of second order that have polynomial eigenfunc- 
tions 1^ (see also |T^). 



The quasiclassical version of Bethe equations (^) describes roots of the 
polynomial eigenf unctions and algebraized part of the spectrum: 

g.(.)-(,-i/2)g;w _, ^ 

E = AhiY. ^nZm + {hf + hf) £ + const. (67) 

n<m m=l 

Diagonalization of (|65|) in terms of Bethe equations has been proposed in 
ref.g. 

If (54(2;) has 4 simple roots, the spectral problem for ( |65D is equivalent to 
the Heun equation p4| (the algebraic forms of Lame and Mathieu equations 
are degenerate cases of the Heun equation). We, therefore, may call the 
difference equation (|65D as g-Heun (Mathieu, Lame) algebraic equations. 

There is also degenerate (triangular) case. The case (ii) corresponds to 
the hypergeometric differential operator 

Q2(^)^ + - (j - l/2)g'2W)^ • (68) 

Its eigenfunctions are classical orthogonal polynomials. The classical limit of 
operators of the case (iii) of Sect. 3 are first order differential operators. 

At last, let us discuss the continuum limit of the linear form (|20|) . Adopt- 
ing the similar notation for ^-dependent coefficients (a = a*-"-* + UaP"^ , etc) as 
^ and putting a*-"-* = S^^ b^^^ = c*-"-* = we obtain the operator: 

L = 2h\2a^°'>Sl + PiSi) , (69) 

where P3 = a^^^ — (3+ = b^^\ (3- = c^^\ As a differential operator it is: 

(P d 

z''-^ + (-/u^' + (1 - 2j + u)z + A)— + 2J^iz, 

where ^i = feW/(2a(°^), = (^^'^ - d^^^) / i2a^'^^) , X = c«/(2a^°^). After the 
change of variable z = and a "gauge" transformation eliminating the 
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term with first derivative, this operator aquires the form ^ — V{x) with 

V{x) = \{fi'e-'^ + X'e'n - i(/x(2j + 1 + u)e-^ + A(2j + 1 - z.)e^) . 
If c^^^ = or b^^^ = it is the Morse potential. 

6 Difference equations solvable in symmetric 
Laurent and Askey- Wilson polynomials 

In this section we describe a related family of difference operators whith 
eigenfunctions given by symmetric Laurent polynomials in one variable y, 
i.e., invariant under the inversion y y^^. Differential operators of this kind 
may be obtained from (^) and (|68|) by the change of variable z = {y + y~^)/2 
(that is the change z = cos 6, e*^ = y which turns Legendre polynomials into 
spherical harmonics). In the case of difference equations, there is no direct 
g-analog of " changes of variables" . We must use another operator algebra. 

The proper operator algebra has been obtained in the paper ||35| as a 
certain degenerate case of the Sklyanin algebra pD|: 



DC = qCD, CA = qAC, 
AB - qBA = qDB - BD = -(1/4) (g^ _ q-^){DC - CA) 
[A,D] = {l/A)iq-q-yC 

q-q 1 



-1X3^2 (70) 
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It can be realized by the shift operators (p!8|): 

A=^Ll_{yT^-y-'T_), D = ^^^{-y-'T^ + yT^ 

y-y y-y 
B = —A-{<f^'v' + iW-'y-^ + i)T++ 



(71) 



+(g-'^y-' + g)(g'^-y + i)T-), 

2 

The representation space of this algebra is spanned by y'' + y~'', k =0, l,...,2j 
(symmetric Laurent polynomials). The standard quantum algebra Uqi^sh) 
can be obtained as a contraction of (7C): C = e'^C, A = eA, D = eD, 
B = B, e 0. There are two central elements: 

no = AD + Uq-q''rc', (72) 

^^ = '-^^^ + BC + li^ + ^-')C'- (73) 

In the continuum limit one obtains U{sl2) in the representaion ([T9|) , where 
z = {y + y^^)/2. Homogeneous bilinear (and linear) forms in A, B, C, D give 
rise to a class of difference equations partially solvable in symmetric Laurent 
polynomials. 

The linear form (pO|) (with A, B, C, D in place of A, B, C, D) gives the 
operator 

2{y-y-^) + 2{y-y-^) ^ > 

where 

a{y) = -b{y^ + y'^) + 2ay - 2dy~^ + 4c - h{q + q~^). (75) 
If 6 7^ and c 7^ the algebraic eigenfunctions of ( [74|) are given by 

'^{y) = y-^^\{{y~yi){y-yT') (76) 
1=1 
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and y/'s satisfy the Bethe equations 

a{q-^yk){qyk - Vk^) i=i^i^k ivk - m) ivk - qyi^) ' 

with the eigenvalues being expressed through y^'s as follows: 

E = aq^ + dq~^ - ^6(g - q-') J^iVi + Vi') (78) 



1=1 



(81) 



(compare with (^51) . 

The general bilinear form ( pUj ) in A, B, C, D gives the difference operator 

G = A{y){Tl - 1) + A{y-'){Tl - 1) + W{y), (79) 

where 

[y-y ^){qy-q ^y ^) 
Wiy) = liq'^ - q-'^) (6i(g2^-i - q-'^+')iy' + y-')+ 

+2{h2q' + hq~^){y + y-^)), 

and 

ao = 4ci + \hi{q^ + g-^ + i)^ = 2g(c2 + K?~^^2 - (g + g"^)fc3)), 
a-i = -2g-Hc3 + - (g + g"')&2)), ^2 = + \q{q + g-')'&i, 
a_2 = q-M + \q-\q + q-^h, 
% = -^g^a, a_3 = \q~^h2, = Whi, a_4 = ig^^&i- 



(82) 



In the generic case (at least one of 6j's is not 0) the algebraic eigenfunc- 
tions of ( |79D are given by ([76|), where yis obey the following system of Bethe 
equations: 

Es=-4 asq'^'yk' ^ X {qhk-yi)iqhkyi - i) 
Et=-4(^sq-^'yt i={j^k hjk - q^yi){ykyi - q^) ' 
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The eigenvalues are 

l<m 

2^' (84) 

- q-'){b2q-'^' + hq^-')Y.{yi + VT^) + {q'' " Q-'')^^ - d)- 

1=1 

-\hi{q + q~'W - q^'^? + \jhi{q - q~'){q^ - q-'). 

The "exactly solvable" case 6i = 62 = ^3 = when the operator is 
triangular is of particular importance. Its eigenfunctions are known as Askey- 
Wilson polynomials p2|, |23|. In this case W{y) = in (|79D . For references 
we recall the conventional form of the difference equation for Askey- Wilson 
polynomials: 

A{y){^>n{q'y) - ^>n{y)) + A{y-'){^>n{q-'y) - = 

(85) 

= (g 2" - 1)(1 - '^WiW2W3W4)'^n{y) 

with 

^iy) = ,F'--triL (86) 

(1 -r)(i - rr) 

where Wi are independent parameters^. The zeros yi, yf^ of the Askey- Wilson 
polynomials satisfy the system of Bethe equations (|83|). In the notation of 
eqs. (p5|) ,(|86D they are 

yr yk-Wg ^ " {q^Vk - yi){q^ykyi - i) 
iiiw^ai/fc-i iJ:}^k (yk- Q^yi){ykyi- q^) 

The form of these equations suggests to ask for an interpretation of the Askey- 
Wilson polynomials in terms of integrable spin chains with boundaries. 

Setting y = e^^ , wi = q"' , W2 = q'^ , W3 = —q'^, W4 = —q^, one readily finds 
that in the continuum limit the difference operator in the l.h.s. of (|85|) turns 
into 

^ :2(7 + 5-a-/3) + (2-a-/3-7-5)cosh(2x))-^. (88) 



dx"^ sinh(2x) dx 



^They are connected with the previous ones in an obvious way; in particular, 
W1W2W3W4 = q^~'^^a/d. 
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Note that in the continuum hmit the number of independent parameters is 
reduced to 2. Ehminating the first derivative term by means of a suitable 
"gauge" transformation, one arrives at the Schrodinger equation for the gen- 
erahzed Poschl- Teller potential. The same potential appears from (|68D after 
a suitable change of variable and a "gauge" transformation. The solutions of 
the Poschl- Teller equation are Jacobi polynomials. 

A particular case of Askey- Wilson polynoms Wi = —w^, W2 = —W4 = g is 
Rogers- Askey- Ismail (or g-Gegenbauer) polynomials p3|. The corresponding 



spectral problem (|85|) in this case have a simple multivariable generalization. 
It is the g-analog (or relativistic extension) of the Calogero-Moser system 

The interrelation between big g- Jacobi and Askey- Wilson polynoms is 
two- fold. First, the former can be obtained p6| under a certain scaling limit 
from the latter, although no parameter is lost in the course of this limit. The 
two families appear to be equivalent. Furthemore, there exists a similarity 
transformation (isospectral transformation) connecting the two families (see 
the end of the next section). One can interpret this transformation as change 
of variables in the difference equation. 



7 Conclusion and discusion 



In this paper we have attempted to classify the difference and discrete op- 
erators that preserve a space of polynomials. These operators may be rep- 
resented as bilinear (or linear) forms of the generators of Uq{sl2) (q-analogs 
of the quantum Euler top). Difference and discrete equations appearing as 
spectral problems for these operators have polynomial solutions. Roots of 
polynomials and eigenvalues are given by the Bethe Ansatz equations. Gener- 
ally polynomial eigenfunctions cover only a part of the spectrum. Sometimes 
all the spectrum is algebraized. Differential equations having polynomial so- 
lutions appear in a certain continuous limit of the difference equations. They 



have been classified previously in Refs. |jTO|, ||, [TT . 

Some periodic difference equations with an incommensurate period have 
a peculiar multifractal spectrum (Harper's equation (|1^) is an example). The 
Bethe Ansatz equations may be helpfull for a statistical description of the 
multifractality. 
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Algebraized difference equations have an intimate relation with integrable 
quantum and classical nonlinear systems. Below we discuss a particular 
aspect of this relation. 

The Heun (Mathieu, Lame) operator (^) written in the algebraic form 
may be transformed to the Schrodinger-type) operator 

H=^,+ V{x) (89) 

by the change of variable 

(90) 



and a subsequent "gauge" transformation. This is the transcendental form oi 
the Heun (Mathieu, Lame) operator. The potential V{x) is a certain elliptic 
function Q In the case of Lame equation it is the Weierstrass function: 
V{x) = n{n + l)p{x). 

Correspondence between solutions of the algebraic and transcendental 
forms of the equations is peculiar. The potential V{x) is a periodic function, 
so the spectrum of the operator (|^) has band structure. Some polynomial 
eigenfunctions of the algebraic form ( 165|) describe edges of the bands, others 



correspond to unnormalizable eigenfunctions |jT0l,llT2[. Both operators (q^) 
and (|65D may have also other eigenfunctions which do not correspond to each 
other. 

An important property of the potential of the Lame equation in the tran- 
scendental form is its direct generalization to the many body case: 

i i>j 



This is the elliptic version of the celebrated Calogero-Moser model . Than 
several questions arise: 

(i) Is there a many-body version of the Lame equation in the algebraic 
form? If yes, it must be a new quantum integrable theory, equivalent to the 
Calogero-Moser model, and having an explicit quantum group symmetry; 

(ii) Do the general Heun (Mathieu) operator in algebraic (|65D or tran- 
scendental form have a many-body generalization? 



^For its explicit form see [D, 
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The same questions apUy to the g-deformed equations: 

(iii) The operator ( pOD is the g-deformation of the algebraic form (|65|) of 
the Heun operator. What it is a g-analog of its transcendental form (jS^) and 
the transformation (pOD? 

(iv) The g-deformation of the Calogero-Moser hamiltonian is known. It 



is the Macdonald operator What is the g-deformed many-body version 
of the general Heun operator? What is the algebraic form of this many-body 
theory? 

Here we may give a partial answer to the question (iii). We suggest that 
the difference equation for the Askey- Wilson polynomials |^ (see eq. 



) is the g-deformation of the trigonometric limit of the Lame operator in 
the transcendental form. Under some a similarity transformation it becomes 
equivalent to the triangular operator of the case (ii) of Sect. 3. The explicit 
form of this transformation is quite complicated. Its existence follows from 
the representation theory of the quadratic algebra (^) developed in the paper 
29[]. We plan to address these questions elsewhere. 
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Appendix A 



Let us show that hnear difference operators of the form (|T]) have (n + 1)- 
dimensional invariant subspace of polynomials Poln iff the coefficients a{z), d{z), v{z) 
are certain Laurent polynomials of order 2 with 7 independent parameters. 
We assume that n > 1. Let us denote fi{z) = a{z), f2{z) = v{z), fsiz) = 
d{z), and consider the decomposition of the coefficients 

m = ft^ + ft\^) + ft\z). 

into a constant fi^\ strictly positive and negative Z/^-* degrees of z. 
Similarly, the operator A = A^°) + A^+) + A*^^-*, where 

A = /iT+ + /2 + fsT^ 

is decomposed into a raising (lowering) parts A'^+)(A'^^) and a diagonal oper- 
ator A^^\ Clearly, A leaves Po/„ invariant if and only if both A'^+) and A^~^ 
do so separately. Consider, say, A(+\ To preserve the space Poln it must 
annihilate and transform each vector z^ {k < n) to a linear combination 
of 

^fc+i^ z^"*"^,...,^". This gives a number of conditions, the first three of them 
are 

q-~^fi-^\z) + ft\z) + q'-^'ft\z) = A^^z, (Al) 
q-'f[''\z) + ft\z) + q-^'ft\z) = A^^^z^ + A^h , 

where A^j^^ are arbitrary coefficients. This is a system of linear equations for 
fi'^^ with non-zero determinant (it is equal to {q — q~^){q + q~^ ~2)). There- 
fore, fi~^\z) are uniquely determined from (|A1|) to be linear combinations 
of z and with 3 independent parameters a\'^\ All the other conditions 
are then automatically satisfied. Similar arguments applyed to A^~^ allow 
one to find //"^ {z) as linear combinations of 2; ^ and z ^ parametrized by 3 
independent constants a[ \ 

The total number of independent parameters is 9 {A\'^\ f^^^; i =1,2,3), 
i.e., the linear space of 2-nd order difference operators preserving the space 
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of polynomials is 9-dimensional. Two of these parameters correspond to 
the constant term and the common factor in ([T8|), so there are 7 essential 
parameters (not including q). 

Calculating fl^\z) explicitly, one finds that A ( [TE| ) coincides (after ob- 
vious changes n —>■ 2j , q ^ q^) with the operator (|30|) obtained by means of 



Appendix B 

Bl. There are other weight representations of Uq{sl2) by the "shift" operators 



T± (see (|T8|)). One of them is 

A = q-^T+, D = q^T_, 

B = {q- q-^)-^z{-q-^T+ + q^ - q-^-^ + q^-^T_), (Bl) 

C = {q- q-^)-^z-\q-^T+ + q{q^ - q-^-^) - q^+^TJ) . 

The representation space again consists of polynomials of degree 2j and the 
value of the Casimir operator is the same as in ([T7|) . The classical limit of 
( [BT| ) is the same as that of (|T^). Bearing in mind applications to 2-nd order 
difference equations, ([191) is more convenient than (|B1|) because any homo- 



geneous bilinear form in A, B, C, D realized as in (|T9D becomes a 2-nd order 
difference operator, whereas a general bilinear form in the representation 
^[l]) gives an operator of fourth order. 

Another representation is realized in the space of Taylor series in z PP| , ^]: 



A = q^T., D = q-^T+, 

B={q- q-')-^z-\q-^^-\l - Tl) + g2j+i(l _ T^)), (B2) 
C = z. 

There is an invariant subspace spanned by z^ with k > 2j. The induced 
representation in the (finite-dimensional) factorspace is equivalent to the spin 
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j representation. The representation ( |B2| ) (in a slightly different form) was 
used in [39| for discretizing the Schrodinger equation with x~^-potential. 



B2. Algebraized difference equations may be obtained from any weight 
representations of g-deformations of s/2 other than Uq{sl2) (p^. 



One of them is generated by J±, Jq EH- 

q'^J+Jo - qJoJ+ = -J+, qJ-Jo - q'^JoJ- = J- 

q-^J+J- - q^J-J+ = iq + q-^)Jo- 
This algebra covers a part of Ug{sl2): 

= BD = [q- q-^)-^qh{q^^Tl - q-^^), 

j_ = CD = {q- q-^)-\^z-\l - Tl), 



(B3) 



(B4) 



Jo = iq-q-T\^-wD 



q + q 



1 



where 



w 



{q + q'')-\q - q-')'^ 



Another "half of the algebra may be obtained by the authomorphism: A 
D ^ A, B ^ C ^ B. 



This algebra has been used in the papers |]T3|, |I0|, ^ to construct differ- 
ence equations having polinomial solutions. A class of difference equations 
obtained from this algebra is of the form = eT^'^, where G is the in- 
vertable Jacobi operator (^91) . Generally these equations do not correspond 
to any hermitian operator. 

All other deformations of s/2 have a similar embedding into Uq{sl2). More 
information about different g-deformed algebras and interrelations between 
them may be found in |Q. Some particular examples were studied in 



[ pm , p2| . Another algebra of difference operators which preserves the space 
of symmetric Laurent polynomials is considered in Sect. 6. 
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